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We study the capture of galactic dark matter by the Solar System. The effect is due to 
the gravitational three-body interaction between the Sun, one of the planets, and a dark 
matter particle. The analytical estimate for the capture cross-section is derived and the 
upper and lower bounds for the total mass of the captured dark matter particles are 
found. The estimates for their density are less reliable. The most optimistic of them give 
an enhancement of dark matter density by about three orders of magnitudes compared 
to its value in our Galaxy. However, even this optimistic value remains below the best 
present observational upper limits by about two orders of magnitude. 
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1. Introduction 

The Dark Matter (DM) density in our Galaxy is about (see e.g. Ref. [T]): 

~ 4 • 10-25 _ 

However, only upper limits on the level of 10^^^ g/cm^ (see below) are known for 
the density of Dark Matter Particles (DMP) in the Solar System (SS). Meanwhile, 
information on their density is of great importance for the experiments aimed at 
the detection of DM. 

The question of capture of Weakly Interacting Massive Particles (WIMP) or 
DMP by the SS was pioneered in Refs.[2j[3l Very recently these studies are pushed 
further by extensive numerical simulations performed in Refs. 131 [5l Our interest to 
this problem was attracted by a recent Ref. £, where an estimate is given for the 
DM density in SS, as resulting from the gravitational capture of galactic DMP's. 
According to the conclusions of Ref. 6j the density of the captured DM, for instance. 
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at the Earth orbit is about 10^^" g/cni'^, only about an order of magnitude below 
the best upper limits on it. 

In the present note we perform new analytical analysis of the gravitational cap- 
ture of galactic DMP's by the SS. According to our results, the increase of the DM 
density in the SS due to this capture is small, certainly well below 10~^° g/cm^. 

2. Dimensional estimate for mass of captured dark matter 

The Solar System is immersed in the halo of dark matter and moves together with it 
around the center of our Galaxy. To simplify the estimates, we assume that the Sun 
is at rest with respect to the halo. The dark matter particles in the halo are assumed 
to have in the reference frame, co-moving with the halo, the Maxwell distribution 
(see Ref-lll): 



with the local rms velocity u ~ 220 km/s. 

Let us elucidate what looks to be the most efficient mechanism of the DMP 
capture. It was pointed out and partly analyzed (though for the capture of comets, 
but not of DMP's), by Petrosky in Ref. 8 and Chirikov and Vecheslavov in Ref. [HI 
Of course, a particle cannot be captured by the Sun alone. The interaction with a 
planet is necessary for it, this is essentially a three-body problem of the Sun, planet 
and DMP. Obviously, the capture is dominated by the particles with orbits close 
to parabolic ones with respect to the Sun, and with the distances between their 
perihelia and the Sun comparable with the radius of planet orbit rp. 

The capture can be effectively described by the so-called restricted three-body 
problem (see for instance Refs. [lOl [TT|) . In this approach the interaction between 
two heavy bodies (the Sun and a planet in our case) is treated exactly. As exactly is 
treated the motion of the third, light body (a DMP in our case) in the gravitational 
field of the two heavy ones. One neglects however the back reaction of a light 
particle upon the motion of the two heavy bodies. Obviously, this approximation is 
fully legitimate for our purpose. Still, the restricted three-body problem is rather 
complicated, and requires in the present case both subtle analytical treatment and 
serious numerical calculations (see Ref. '8) . Under certain conditions the dynamics 
of light particle (e.g. DPM) becomes chaotic. 

However, the amount of DM captured by the SS can be found by means of 
simple estimates. The total mass captured by the Sun (its mass is M) together with 
a planet with mass nip, during the lifetime 





T ~ 4.5 • 10^ years ~ lO" s 



(3) 



of the SS, can be written as follows: 



ArUp = pgT < av > ; 



(4) 
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here a is the capture cross-section. The product av is averaged over distribution ^ ; 
with aU typical velocities in the SS much smaller than u, this distribution simplifies 
to 



To estimate the average value < av >, we resort to dimensional arguments, 
supplemented by two rather obvious physical requirements: the masses nip and M 
of the two heavy components of our restricted three-body problem should enter the 
result symmetrically, and the mass nid of the light component (DMP) should not 
enter the result at all. Thus, we arrive at 



here k is the Newton gravitation constant; an extra power of tt, inserted into this 
expression, is perhaps inherent in a. The final estimate for the captured mass is 



Since the capture would be impossible if the planet were not bound to the Sun, it is 
only natural that the result is proportional to the corresponding effective " coupling 
constant" knipM. 

Thus obtained values for the masses of DM captured due to the planets of the 
SS, are presented in Table 1. We quote also therein the corresponding results of 
Ref. [6]for these masses. The disagreement is huge for all planets, and especially for 
the light ones where it exceeds two orders of magnitude. We cannot spot exactly its 
origin since the calculations of Ref. (61 involve rather complex numerical simulations 
(it is possible that their assumption of capture radius rt, ^ rpiinp/M)^/'^ does not 
correspond to reality). On the other hand, however, we cannot see any reasonable 
possibility for a serious increase of our results. Moreover, in a sense they can be 
considered as upper limits for the amount of the captured DM, at least because 
we have neglected here the inverse process, that of the ejection of a captured DMP 
due to the same three-body gravitational interaction. The result ([7|) is given for the 
three body problem. The dynamical mechanism of capture is described below in 
next Section. The contribution of the diffusive (non-dark) matter in the SS should 
be significantly smaller since a homogeneously distributed dust gives compensation 
of gravitational forces acting on DMP. 

The total mass Amy of the DM captured by the planets is strongly dominated 
by the heavy Jovian planets, Jupiter, Saturn, Uranus, Neptune, and constitutes 
according to Table 1 about Amr ~ 1-5 • 10^^ g. This value is small as compared 
to the total mass ~ 10'^^ g of the common matter in the SS. It is small even as 
compared to the total non-captured mass of the DM in the SS: this total mass, 
calculated with value ([1]) for the DM density, constitutes ^ 10^^ g (we assume 
here that the effective radius of the SS is about 10^ an). However, it is an order 




(5) 




(6) 




(7) 
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of magnitude larger than the DM mass of density ([T|) inside the radius of Neptune 
orbit rN w 30 au. 

The contribution to the discussed effect of the diffuse (non-dark) matter in the 
SS should be significantly smaller since in a homogeneous dust the gravitational 
forces acting on DMP are compensated. 

The dynamical mechanism of capture is described in the next Section. 



Table 1. DM mass captured by planets (in g) 



Planet 


Mercury 


Venus 


Earth 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


this work 


0.22-10^* 


3.2101** 


3.9-10^® 


0.42-10^8 


1239-10^** 


372-10^** 


57-10^* 


67-101* 


Xu, Siegel Ref. |6] 


0.42-1020 


3.5-1020 


3.8-10^0 


1.2-10^0 


49-10^0 


28-10^0 


12-10^0 


16-10^0 



3. Dynamical approach 

For the restricted three-body problem, in a close similarity to dynamics of comets 
(see Refs.lHllll, the DMP dynamics can be described by a symplectic area-preserving 
map 

w = w + F{(j)) , 4>^(j) + 27ru;-3/2 . (8) 

Here w = —2Erp/kmdM is the energy of DMP with mass rud rescaled by its 
gravitational energy on distance rp from the Sun; (/) is the phase of the planet on 
its circular orbit at the moment when the DMP is at the perihelion and F{(j)) is 
a certain periodic function of (j)- Bars denote the new values of variables after one 
rotation around the Sun. The physical meaning of this dynamical map is rather 
simple: the first equation gives the change of DMP energy after one passage near 
the Sun, the second equation gives the change of the planetary phase between two 
passages of DMP and is essentially determined by the Kepler law. The first equation 
is valid also for scattering particles with positive energy {w < 0). Thus DMP can 
be captured by the Sun and the planet only if its rescaled energy \w\ < F^^ax- 
After the capture, the DMP dynamics is described by map ^ until ejection. To 
compute the captured DM mass Am^ we assume that, after being once captured, 
the DMP remains captured for the whole life time T of SS. In this way we obtain 
the maximum bound for Awp. 

The kick function F(0) was computed in Ref. [8] for the case when a comet 
(or DMP) and a planet move in one plane and when the perihelion distance 
q > rp. In this case F{(t)) — {mp/M)(3{q/rp) sin (j) and the function /3(x) w 
26exp(-4a;3/2/3%/2)/2:i/'* so that /3(1) w 10. Effectively the function F is deter- 
mined by the frequency Fourier component of the force between the planet and 
DMP, since the rotation of the planet is rapid compared to the rotation of DMP 
the amplitude of the component is exponentially small for q ^ rp when the DMP 
motion is smooth and analytical. In this case /3 is exponentially small and there is 
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practically no trapping of DMP. For q ^ Vp the motion is not analytic due to close 
passage between the planet and DMP and /3 is relatively large. In this case a DMP 
with rescaled energies —w < Pirip/M can be captured by the planet. It is interesting 
to note that the map ([8]) with F{(j)) ~ sin is known as the Kepler map. It describes 
the process of microwave ionization of Rydberg atoms and chaotic autoionization 
of molecular Rydberg states (see Ref. [T2|and Refs. therein). 

We note that the energy change of DMP given by F{(j)) results from the inte- 
gration over the whole orbit rotation of DMP around the Sun which includes many 
orbital periods of the planet. Thus this energy change appears from long-range in- 
teraction and has qualitatively different origin compared to local close collisions 
between DMP and planet which were assumed to give the main contribution for 
DMP energy change in Refs. [21 [3l [6] 

Let us now estimate the capture cross-section a assuming that for all DMP the 
dynamics is described by the Kepler map with fixed (3^1, Then only DMP with 
energies \w\ = v^Tp/knipM = v^/Wp < Pnip/M are captured under the condition 
that q < Tp (here Vp is the velocity of the planet). The value of q can be expressed via 
the DMP parameters at infinity where its velocity is v and its impact parameter is 
Td and hence q = (vrd)'^ /2kM (see Ref. [T3| . Since q ^ rp we obtain the cross-section 

2TTkMrp/v^ ~ 2TTrl{vp/vf - 2TTrlM/{pmp) , (9) 

where the last relation is taken for those typical velocities ~ /Svpirip/M at which 
the capture of DMP takes place (for q « lAvp we have f3 « 5). Then Eqs. (g]), ^ 
give the captured mass Amp of ([7]) with an additional numerical factor /3 ~ 1. 

According the above estimates DPM captured by Jupiter have typical velocities 
at infinity v ~ {[3mp/MY^'^Vp ^ 1km/ s for typical /3 ~ 5 and mp/M « 10~^, 
Vp « 13/cm/s. This value of w is in a good agreement with the numerical simulations 
of Peter Ref. [5] which give typical captured DPM velocities for Jupiter of 1km/ s. 



Table 2. Density of DM for planets (in g/cm^) 



Planet 


Mercury 


Venus 


Earth 


Mars 


Jupiter 


Saturn 




2.7-10-22 


6.0-10-22 


2.7-10-22 


8.4-10-2-1 


6.2-10-22 


3.0-10-23 


Pdm 


1.810-21 


1.5-10-21 


9.3-10-22 


6.6-10-22 


6.5-10-22 


3.0-10-23 



Table 3. Angle of perihelia precession (in seconds per century) 



Planet 


Venus 


Earth 


Mars 




8.6248 


3.8388 


1.3510 


&<l>ohs 


8.6247 ± 0.0002 


3.8390 ± 0.0003 


1.3512 ± 0.0003 



Another interesting feature of the analytical expression for the cross-section 
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of captured particles u (H)) is that it is much larger than the area of the planet 
orbit. In fact, a diverges at small velocities as a ~ but this divergence is 

weaker compared to the divergence of the Reserford cross-section. In our case of 
the restricted three-body problem the divergence appears due to the property of 
the Kepler motion where the DMP distance at perihelion is proportional to square 
of the orbital momentum which in its turn is proportional to the product of the 
velocity v and impact parameter at infinity. In addition it is important to use the 
value of typical DMP velocity captured by the planet for perihelion distance of the 
order of rp. This leads to the analytical equation ([9]) for the capture cross-section 
in the restricted three-body problem. 

4. Density of dark matter 

While the total masses Amp of the captured DM can be (hopefully) described by 
simple dimensional estimate (l7|), the situation for the corresponding DM densities 
App is more subtle. The reason is as follows. The captured DMP's had initial trajec- 
tories predominantly close to parabolas with respect to the Sun, and the velocities 
of these DMP's change only slightly as a result of scattering. Therefore, it is quite 
natural that the final, elliptical trajectories of these DMP's have large semi-major 
axes. 

Indeed, DMP captured into an elliptic trajectory had initially a hyperbolic tra- 
jectory, focussed at the Sun and close to a parabolic one. As a result of the capture, 
the eccentricity e of the trajectory changes from e = 1 + ei to e = 1 — 62, with 
e = ei + 62 <C 1. It is quite natural that the final, elliptical trajectories of the 
captured DMP's have large semi-major axes. 

To estimate their typical values, we recall (see Ref. [13]) that the radius-vector 
r of a captured DMP (counted off the Sun) is related to the azimuthal angle 4> as 
follows: 

r = p/(l + ecos(/)) , (10) 

where p is the so-called orbit parameter (its value is irrelevant for our line of rea- 
soning). Obviously, the maximal rmax and minimal distance rmin from the Sun 
correspond to cos0 = ±1, so that their ratio is 

rmax/rrnin = {I + e)/{l- e) . (11) 

In the numerator of this ratio, we can safely put with our accuracy 1 -I- e ~ 2. As 
to its denominator, we recall that the difference 1 — e is related to the gravitational 
perturbation by planet, and therefore is proportional to nip. Thus, by dimensional 
reasons, 

max I min ^ Ad/uip . (l'^) 

The minimal distance between DMP and the Sun rmin should be on the same 
order of magnitude as the radius rp of the planet orbit. Therefore, the semi-major 
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axis a of resulting ellipse is huge: 

(M/nip) . (13) 

In particular, in the case of Jupiter our estimate gives ^ lO^rp. A similar numer- 
ical factor appears in Ref. [51 

According to the numerical calculations of Ref. [8l in the case of Jupiter the 
values of the semi-major axes for the resulting trajectories belong to the interval 
10'^ — 10^ au for q/rp = 4 — 6. The fact that it is comets that are considered 
in Ref. [51 but not DMP's, is obviously of no importance for this conclusion. The 
minimum value of ra is defined by the maximum Wch ^ Tpjra value which can 
be reach by an injected DMP during its chaotic motion. In fact Wch is the chaos 
border and according to the Chirikov resonance-overlap criterion (see Ref. [M]) we 
have Wch ~ {i'K(3{mp/MYf'^^ as it was shown in Ref. 8. For Jupiter nip/M 10"'^ 
and at /3 « 5 corresponding to q « 1.4rp we have Vajvp « \/wch ~ 3. We note that 
this value of (3 gives the maximum Fmax ~ 0.005 corresponding to the similar value 
foimd for the comet Halley (see Ref. [HI. In fact, the data presented in Ref. 9 show 
that the comet Halley have the chaos border around Wch ~ ^p/va ~ 0.3 (see Fig. 3 
in Ref.H. 

Of course the values of ra linked to the chaos border in w are the minimum 
ones since during its chaotic dynamics DMP have also < ui ^ Wch with larger ra ■ 
However, we are interested in orbits captured for very large times T Such times 
are by two orders of magnitude larger than a typical diffusive life time of comet 
Halley found to be of the order of 10^ years (see Ref. [9]). It is known that chaotic 
trajectories may be sticking to boundaries of integrable islands for very long times 
(see Ref. [15] and Refs. therein) and hence we can expect that those orbits will be 
somewhere in vicinity of the chaos border around Wch ~ 1 with ra ^ rp. In fact, for 
the case when the inclination angle between the planes of DMP and planet 9i > Q 
and when q < rp the function i^(0) contains higher harmonics of (see the case of 
the comet Halley in Ref. [HI). This leads to easier emergence of chaos so that even 
for light planets one may have the chaos border Wch ^ 1- 

Therefore we can make an assumption resulting in the most optimistic prediction 
for the "partial" dark matter densities Apj, . We assume that each of the total masses 
Amp of the captured DM occupies the volume (47r/3)rp where rp is the orbit radius 
of the corresponding planet. We do not claim that this assumption is correct, but 
believe, however, that the comparison of its (almost certainly, overoptimistic) results 
with the observational limits will be instructive. The corresponding values of the 
"partial" dark matter densities App — Amp/{4:TTrp/3) (in g/cm^) are presented 
in Table 2. We omit in it the densities due to Uranus and Neptune, tiny even at 
the discussed scale. Then, in accordance with the accepted model, the total dark 
matter density puM at a given radius does not coincide with the corresponding App. 
It includes, in line with it, the sum of the contributions to the density due to all 
the planets, outer with respect to the given one. 
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5. Ergodic time scale 

The estimates given above neglect ejection of DMP from the SS. Such an assumption 
is not justified if the DMP dynamics in SS becomes completely ergodic on a time 
scale Te <C T. Then after the time Te the detailed balance principle becomes valid 
and the density of captured DMP becomes as its galactic density as it was argued in 
Ref. [2|(see also discussion in Refs.[3j|4l[5]). However, the estimate of Te have certain 
subtle points. In the frame of the map ([8]) it is given by the diffusion time from w = 
to w — Wch- For the Kepler map the diffusion coefficient is D « /3^(mp/Af)^/2 and 
hence ~ 2{Mwchl PrUpfTd where Tj, is an average period of DMP. For the case 
of Jupiter such an estimate gives a satisfactory value of Te ~ 10^ years for the 
case of the comet Halley as it is discussed in Ref. [H However, /3 sharply decreases 
with the increase of perihelion distance g oc where £ is the orbital momentum 
of DMP. As a result a growth of ^ can give sharp increase of which can become 
comparable with T. The effects linked to variations of i were not considered in 
Ref. ini They are properly treated in the numerical simulations of Refs. H] [51 but 
there only the effect of Jupiter is considered. Other planets and fluctuating galactic 
fields can give stronger growth of £ with a signiflcant increase of T^. Indeed, from the 
studies of Rydberg atoms in a microwave field it is known that time-oscillating space 
homogeneous fields can produce strong variation of eccentricity of orbits (see Ref. [TBI 
and Fig. 11 there). Also it is known that noise generates penetration of chaotic 
trajectories inside integrable islands and very slow decay of Poincare recurrences 
with diverging trapping time (see Ref. I17p . In addition to that recent large scale 
numerical simulations of Ref. [T8| show significant changes of eccentricity of planets 
on a time scale of the order of T. Therefore the question about value for DMP 
captured by the SS requires further studies. It is not excluded that it is comparable 
or even larger than T . In such a situation the upper bound ([7|) will be close to the 
real value of total captured mass. 

In any case it is clear that there is practically no ejections of captured particles 
on a time scale of DMP orbital period Tc- A typical captured DMP rescaled energy 
is Wc ~ (3mp/M corresponding to one iteration of the map ([8|) which gives a change 
of w from negative to positive values. The rotation period of such DMP is rather 
large compared to the period Tp of the planet: Tc/Tp ~ {j3mp/M)~'^/'^ . For the case 
of Jupiter Tp « 11 years and at /3 ~ 1 we have Tc ~ 3 • lO'^Tp ^ 3 • 10^ years. This is 
by a factor 10'* shorter than the SS life time T. This gives the lower bound of the 
captured DM mass which is obtained by replacing T by Tc in the equation ([7|). 

6. Observational upper limit on the density of dark matter 

At last, let us consider the observational data on the DM in our SS. The most reliable 
and accurate information on it follows from the studies of the perihelion precession 
of Venus, Earth, and Mars. Under the assumptions that the DM density pdm is 
distributed spherically symmetric with respect to the Sun and that the eccentricity 
of the planetary orbit is small, the corresponding relative shift of the perihelion per 
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period is (see Ref. [TOl and Refs. therein): 



M 



(14) 



where r is the radius of the orbit. This relation gets almost obvious (up to an overall 
numerical factor) if one recalls that, in virtue of the Gauss theorem, for a spherically 
symmetric density p{r) the action of the DM inside the orbit reduces to that of a 
point-like mass, and therefore does not induce the perihelion precession. On the 
other hand, for such density p{r), the DM outside the orbit does not influence at 
all the motion of a planet. 

The recent, most precise observational data (see Ref. i20j) on the precession of 
perihelia are presented in Table 3 (therein the theoretical values S(j)th of the perihe- 
lion rotation and the results of observations (5(/)obs are given in angular seconds per 
century). With these data, one arrives at the upper limits on the DM density at 
the distances from the Sun, corresponding to the orbit radii of Venus, Earth, and 
Mars, on the level of 



This observational upper limit exceeds by about two orders of magnitude the 
results (almost certainly overestimated) presented in Table 2. 

7. Summary 

Our results do not mean, however, that the searches for the dark matter in the SS 
are senseless. Of course, the capture of the Galactic DM analyzed here is not the 
only conceivable source of the dark matter in the Solar System. It is quite possible 
in particular that the Solar System itself has arisen due to a local high-density 
fluctuation of the dark matter. 

Now on the related theoretical problems. To obtain more firm results for the 
captured DM mass and density one needs to take into account that the kick function 
F(0) in the map ([8]) depends on an inclination angle between planes of DMP and 
planet orbits and also on the DMP perihelion distance. However, a typical case of 
the comet Halley analyzed by Chirikov and Vecheslavov in Ref. [HI gives the map 
function of a form similar to one discussed here so that the estimates presented 
should be also applicable for such more general DMP orbits. Further analytical and 
numerical studies are required for a better understanding of DMP dynamics inside 
the Solar System. 
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